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ABSTRACT 

We present the discovery of fifteen new T2.5-T7.5 dwarfs (with estimated distances 
between ~24-93pc), identified in the first three main data releases of the UKIRT 
Infrared Deep Sky Survey. This brings the total number of T dwarfs discovered in the 
Large Area Survey (to date) to 28. These discoveries are confirmed by near infrared 
spectroscopy, from which we derive spectral types on the unified scheme of Burgasser 
et al. (2006). Seven of the new T dwarfs have spectral types of T2.5-T4.5, five have 
spectral types of T5-T5.5, one is a T6.5p, and two are T7-7.5. We assess spectral 
morphology and colours to identify T dwarfs in our sample that may have non-typical 
physical properties (by comparison to solar neighbourhood populations), and find that 
one of these new T dwarfs may be metal poor, three may have low surface gravity, 
and one may have high surface gravity. The colours of the full sample of LAS T dwarfs 
show a possible trend to bluer Y — J with decreasing effective temperature, and some 
interesting colour changes in J — H and z — J (deserving further investigation) beyond 
T8. The LAS T dwarf sample from the first and second main data releases show good 
evidence for a consistent level of completion to J=19. By accounting for the main 
sources of incompleteness (selection, follow-up and spatial) as well as the effects of 
unresolved binarity and Malmquist bias, we estimate that there are 17±4 ^T4 dwarfs 
in the J ^19 volume of the LAS second data release. Comparing this to theoretical 
predictions is most consistent with a sub-stellar mass function exponent a between 
-1.0 and 0. This is consistent with the latest 2MASS/SDSS constraint (which is based 
on lower number statistics), and is significantly lower than the a ~ 1.0 suggested by 
L dwarf field populations, possibly a result of the lower mass range probed by the T 
dwarf class. 
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1 INTRODUCTION 

The advent of the large scale Sloan Digital Sky Survey 
(SDSS; York et al. 2000), the 2-Micron All Sky Survey 
(2MASS; Skrutskie et al. 2006) and Deep Near Infrared Sur- 
vey of the Southern Sky (DENIS; Epchtein et al. 1997) were 

the main factors that led to the identification and study of 
the two spectral classes beyond M; the L dwarfs (with effec- 
tive temperature T^s ~2300-1450K; Gohmowski et al. 2004) 
have dusty atmospheres and very red near infrared colours, 
while the even cooler T dwarfs (Tcff <1450K) have clear at- 
mospheres (from which the dust has settled), and their blue 
infrared spectra are dominated by strong CH4 and H2O ab- 
sorption bands. The spectral typing scheme of Kirkpatrick 
et al. (1999) established the first properly classified L dwarf 
population via 2MASS discoveries. The first T dwarf was dis- 
covered as a close companion to the early M dwarf GL229 
(Nakajima et al. 1995), and the first L dwarf as a compan- 
ion to the white dwarf GD 165 (Becklin & Zuckerman 1988; 
Kirkpatrick ot al. 1993). However, the subsequent growth of 
the known T (and L) dwarf populations (>500 L dwarfs and 
>100 T dwarfs are known at time of writing; Kirkpatrick 
2005) has been predominantly achieved via searches of the 
SDSS (e.g. Leggett et al. 2000; Geballe et al. 2002; Knapp 
et al. 2004; Chiu et al. 2006) and 2MASS (e.g. Kirkpatrick 
et al. 2000; Cruz et al. 2007; Burgasser et al. 2002; 2004; 
Tinney et al. 2005; Looper Kirkpatrick & Burgasser 2007) 
databases. 

The coolest T^a probed by SDSS and 2MASS are cur- 
rently defined by the eight T7.5-8 dwarfs (as typed by the 
Burgasser et al. (2006b; B06) scheme) discovered in these 
surveys. These objects have T^s in the range 725-950K 
(Geballe ct al. 2001; Saumon et al. 2006; 2007; Leggett et 
al. 2007; Burgasser, Burrows & Kirkpatrick 2006). At lower 
Tcff it is possible that a new spectral class (pre-emptively 
called Y, following Kirkpatrick et al. (1999)) may be nec- 
essary, if for example, spectral changes occur such as the 
strengthening of ammonia absorption, and/or the conden- 
sation of atmospheric water clouds (Burrows, Sudarsky & 
Lunine (2003)). L and T dwarfs already encompass the tem- 
perature scale of transiting "hot Jupiters" (e.g. Knutson et 
al. 2007), and populating the even cooler Toef regime will 
allow us to study and understand atmospheres whose Tcff is 
similar to cooler extra-solar giant planets populations. 

The current complement of known L and T dwarfs al- 
lows some constraints to be placed on the substcllar mass 
function (e.g. Chabrier 2003 and references therein; Allen et 
al. 2005; Metchev et al. 2007). However, if the brown dwarf 
mass function is to be accurately constrained, a significantly 
larger number of late T dwarfs will be extremely beneficial, 
since the T^s distribution of T dwarfs in the <950K range is 
particularly sensitive to mass function variations (e.g. fig 5 
of Burgasser 2004) . Also, accurate constraints on the brown 
dwarf birth rate (i.e. their formation history) need the im- 
proved statistics that come with larger numbers of both L 
and T dwarfs. 

The UKIRT (UK Infrared Telescope) Infrared Deep Sky 
Survey (UKIDSS; Lawrence et al. 2007) is a new infrared 
survey being conducted with the UKIRT Wide Field Cam- 
era (WFCAM; Casah et al. 2007). UKIDSS is a set of five 
sub-surveys, with three wide field surveys - the Large Area 
Survey (LAS), the Galactic Cluster Survey (GCS) and the 



Galactic Plane Survey (GPS) - and two very deep surveys - 
the Deep Extra-galactic Survey (DXS) and the Ultra-Deep 
Survey (UDS). The LAS is the largest of the wide- field sur- 
veys, and will cover 4000 sq dcgs of sky in four filter bands, 
going several magnitudes deeper than 2MASS. UKIDSS be- 
gan in May 2005, and at the time of writing there have been 
four ESQ (European Southern Observatory)-wide releases, 
including an Early Data Release (EDR) in February 2006 
(Dye et al. 2006), and three subsequent main data releases: 
Data Release 1 (DRl) in 2006 July covering 190 sq degs 
(Warren et al. 2007b), Data Release 2 (DR2) in 2007 March 
covering 280 sq degs (including DRl; Warren et al. 2007c), 
and Data Release 3 (DR3) in 2007 December, covering 900 
sq dogs (including DRl and DR2). 

These data are providing un-paralleled sensitivity to L 
and T dwarf populations, and our collaboration has begun 
a variety of UKIDSS-based searches for these objects. Here 
wo focus on our search for late T dwarfs and potentially 
new record breaking low- Tcff objects. Previous LAS T dwarfs 
have boon presented by Kendall et al. (2007), Lodieu et al. 
(2007), Warren et al. (2007a) and Chiu et al. (2008), who 
have discovered a total of 13 spectroscopically confirmed 
LAS T dwarfs from the EDR and DRl. This included the 
first T8.5 dwarf ULAS J0034-0052 (Warren et al. 2007a). 

Here we report the discovery of 15 new LAS T dwaxfs 
discovered in DRl, DR2 and some of DR3. Section 2 sum- 
marises our selection criteria for identifying candidate T 
dwarfs as well as possible 400-700K objects. Section 3 de- 
scribes the follow-up photometry we have obtained to iden- 
tify spurious objects amongst these candidates, and Section 
4 presents our spectroscopic follow-up and confirmation of 
the 15 new T dwarfs. Section 5 discusses the spectral mor- 
phology and colour of the new T dwarfs and how these could 
relate to their physical properties, and Section 6 presents 
updated constraints on the size of the LAS T dwarf popular 
tion, compared to theoretical predictions. Section 7 discusses 
some future work and gives our conclusions. 



2 IDENTIFYING T DWARF CANDIDATES 

In this section we describe our photometric search for T 
dwarfs and cooler 700-400K objects (potential Y dwarfs). 
We will focus this description on our searches of DRl and 
DR2, since our follow-up of candidates from these data re- 
leases has a reasonably well constrained level of completion, 
and issues such as contamination amongst the candidates 
will be addressed more usefully. However, the process of T 
dwarf identification is ongoing with DR3, using these same 
techniques. Note also that some DRl results have previously 
been published (Kendall et al. 2007; Lodieu et al. 2007; War- 
ren et al. 2007), but for clarity we here describe our search 
of DRl and DR2 as a whole. 

We base our search on the current knowledge of izJHK 
properties of known objects previously discovered in 2MASS 
and SDSS, as well as on colour trends suggested by the latest 
theoretical atmosphere models for Tea ranges below those 
that are well probed by previous surveys. We also make 
use of LAS magnitudes measured using the new Y filter 
(0.97-1.07/im), which was specifically designed and installed 
in WFCAM to ease the selection and separation of high- 
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redshift quasars and cool brown dwarfs (Warren & Hewett 
2002). 

2MASS T dwarfs have neutral to blue near infraxed 
colours with decreasing T^s (Burgasser et al. 2002), and a 
fairly uniform Y — J (from synthesized colours; Howctt 
et al. 2006). Late T dwarfs are clearly well separated from 
redder L dwarfs by their J — H colour (although there is 
some J — H colour overlap for early T dwarfs), and from 
earlier objects by, in general, both J — H and Y — J colour 
(sec Figure 1). The model predictions suggest that the near 
infrared colours will remain blue {J — H < 0.0, J — K < 0.0) 
for 400-700K. However, while the models uniformly predict 
that the 400-700K dwarf colour sequence should form an 
extension of the known T dwarfs, they differ somewhat on 
their predicted Y — J colour trends. The cloud-free Cond 
models (Allard et al. 2001; Baraffe et al. 2003), the more 
recent Settl models (Allard et al., in preparation) and the 
Marley et al. (2002) models all suggest a blue Y — J colour 
trend, while the Burrows et al. (2003) and Tsuji, Nakajima 
& Yanagisawa (2004) models suggest a red Y — J colour 
trend for this Tes regime. 

In addition, known T dwarfs from SDSS have very red 
{i — z ^2.2) optical colours, and extremely rod optical-to- 
infrared colours {z — J ^2.5; note that here and hereafter 
i and z are AB magnitudes whereas J is on the Vega sys- 
tem), the latter being a particularly good indicator of T^s 
for early-mid T dwarfs (Knapp et al. 2004). The model pre- 
dictions suggest that these extreme optical and optical-to- 
infrared colours should continue to be strong indicators of 
low-Tefj in the 400- TOOK regime. A combination of YJHK 
photometry from the LAS and deep optical constraints from 
SDSS thus offer an extremely powerful tool to identify sam- 
ples of mid-late T dwarfs and even cooler objects to pho- 
tometric depths several magnitudes deeper than previous 
wide-field surveys. 

2.1 Initial sample selection 

Our search methodology was to mine (via the WFCAM 
Science Archive; Hambly et al. 2007) the LAS for objects 
with YJHK photometry consistent with known T dwarfs 
or with the expected colours of 400-700K objects using the 
models as a guide, and then cross-match with SDSS (where 
possible) to obtain optical-to-infrared constraints. Here we 
describe three sets of search criteria that probe to differ- 
ent depths in DRl and DR2, requiring coverage in all four 
YJHK bands, but relying on source detection in different 
numbers of bands. 

Our first search method required YJHK detection with 
y — J ^0.8 and J — H ^0.4 (to separate T dwarfs from L 
dwarfs). This method is limited by the if— band for mid- 
late T dwajfs, does not probe the full range of F — J colour 
for potential 400-700K objects, and is thus better at find- 
ing brighter early T dwarfs which can have redder H — K 
colour. We also cross-matched objects with SDSS, requir- 
ing a detection with i — z ^2.2 and z — J ^2.5. While we 
also performed searches requiring optical non-detection, we 
postpone analysis of these candidates to a future publica- 
tion, and here present the one now object identified with a 
SDSS detection; ULAS J115759.04-F092200.7 (see Table 1). 

Our second search method required YJH LAS detec- 
tion and a K non-detection, and probed the colour space 



shown in Figure 1, including objects with J — H ^0.1 
and Y — J >0.5. This search is _ff-band limited for ob- 
jects with blue J — H, but fully probes the range of 400- 
700K colour space suggested by the models. Candidates were 
cross-matched with SDSS and rejected if found to be optical 
detections with either i — z <2.2 or 2; — J <2.5, or found to 
have no SDSS coverage. This search produced a sample of 
33 candidates. 

Our third search method required YJ LAS detection 
and HK non-detection, with Y — J ^0.5. For objects with 
blue J — H this search is J-band or F-band limited for Y — 
J <0.8 and Y — J >0.8 respectively, and probes to a greater 
depth than method 2. SDSS cross- matching constraints were 
imposed as before, and resulted in a sample of 45 candidates. 

In addition we used our second search method to search 
the LAS sky outside the SDSS DR6 footprint, substituting 
our SDSS optical constraints for shallower Schmidt plate 
/-band coverage from the SupcrCOSMOS sky survey. This 
search identified 24 candidates. 

Figure 1 shows the J — H, Y — J two-colour diagram 
and the J against J — H colour magnitude diagram for the 
T dwarf and Tes =700-400K candidates from DRl and DR2. 
Candidates arc cither crosses (where Y—, ,7— and _ff— band 
detections were available) or arrows (where H was a non 
detection). A sample of L dwarf candidates (plus signs; LAS 
Y - J >0.9, J -H >0.5, J -K >1.2, and SDSS z <20.8) 
is also shown for comparison (some spectroscopically con- 
firmed L dwarfs from this sample will be reported in a sep- 
arate paper) , as is a sample of typical brighter sources from 
the LAS (points). These brighter sources were selected from 
a one degree radius area of DR2 sky at high Galactic lat- 
itude, and have YJHK magnitude uncertainties less than 
0.03 magnitudes and a Sloan counterpart. The L dwarf can- 
didates are being followed up, but we will not discuss them 
further here, showing them purely for comparison. Our YJH 
selection box {Y — J >0.5, J — H < 0.1) is indicated with 
a dashed line, and contains over-plotted boxes (dotted and 
dot-dashed lines) that illustrate the range of model predic- 
tions for reft=700-400K (see figure caption). 

Candidates from our DRl-3 searches that have been 
confirmed as T dwarfs via spectroscopy (sec Section 4) are 
shown in Table 1, which summarises the database release 
in which they were first found, as well as the bands that 
were used to make the searches. All of our new T dwarfs 
are in SDSS DR6 sky, although two previously confirmed 
candidates (Lodieu et al. 2007) arc outside the SDSS DR6 
footprint. Figure 2 presents finding charts for the 15 new T 
dwarfs. 



3 FOLLOW-UP PHOTOMETRY 

Photometric contamination amongst our candidates has sev- 
eral causes. In general, photometry near the survey limits 
will have significant uncertainties due to low signal-to-noise 
ratio and/or survey image defects (see Dye et al. 2006), and 
in certain circumstances normal stars can have their colours 
scattered to the blue. Scattered late M dwarfs with J >18 
will not always be ruled out via SDSS non detection if their 
i - J =3-4 (the 5-0- SDSS i-band limit is ~22). Also, if LAS 
images in different bands are measured on different nights, 
variable sources can yield extremely blue survey colours, and 
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Figure 1. A J — H against Y — J two-colour diagram and a J against J — H colour magnitude diagram showing candidate late T and Y 
dwarfs from DRl and DR2. Y JH(K) candidates are shown as crosses, and Y J candidates as upper limit arrows. An illustrative L dwarf 
candidate sample is shown by plus signs. Our YJH candidate selection box is shown with a dotted line, and contains theoretical model 
predicted colours for Tgf5f=400-700K dwarfs from: the cloud-free Lyon-Cond models (Allard et al. 2001; Baraffe et al. 2003); the more 
recent Lyon-Settl models (Allard et al., in preparation); the AMES models (e.g. Marley et al. 2002; Saumon et al. 2003); the Tucson 
models (e.g. Sharp 2000; Burrows Sudarsky & Hubeny 2006); and the Tsuji models (Tsuji Nakajima & Yanagisawa 2004). 



Table 1. Newly discovered UKIDSS T dwarfs. The table lists the coordinates, the UKIDSS LAS database release in which objects were 
found, and the bands that were used to make the searches in which the objects were identified. All objects have SDSS coverage, although 
only one is detected in SDSS. The table contains the newly confirmed T dwarfs as well as (at the bottom) three that we rule out with 
spectroscopy. 



Name" 


Database 


NIR band 


SDSS band 






detections 


detections 


ULAS .1022200.43 - 002410.5 


DRl 


YJH 




ULAS J082327.46 + 002424.4 


DR3 


YJHK 




ULAS J083756.19 - 004156.0 


DRl 


YJH 




ULAS J085910.69 + 101017.1 


DR3 


YJHK 




ULAS J093829.28 - 001112.6 


DR3 


YJH 




ULAS J093951.04 -1- 001653.6 


DR3 


YJH 




ULAS J095829.86 - 003932.0 


DRl 


YJ 




ULAS J115038.79 + 094942.9 


DR2 


YJH 




ULAS J115759.04 -|- 092200.7 


DR2 


YJHK 


i,z 


ULAS J130303.54 -|- 001627.7 


DRl 


YJ 




ULAS J131508.42 + 082627.4 


DR2 


YJH 




ULAS J150135.33 + 082215.2 


DR2 


YJ 




ULAS J150547.89 + 070316.6 


DR2 


YJ 




ULAS J154427.34 + 081926.6 


DR2 


YJ 




ULAS J154701.84 + 005320.3 


DR2 


YJH 




ULAS J003925.90 + 000257.7 


DRl 


YJH 




ULAS J130150.35 + 002314.8*' 


EDR 


YJH 




ULAS J224355.03 - 004618.5 


DRl 


YJH 





"Here the source names include full coordinates. In subsequent tables/figures the 
source names are truncated. 

''First reported as a T dwarf candidate by Kendall et al. (2007). 
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Figure 2. Finding charts for the 15 new T dwarfs discovered in the UKIDSS LAS and presented in this paper. Charts are J— band 
images of generally two arcmin on a side (unless constrained by the edge of a UKIDSS frame) with north up and east left. Indicated 
spectral types are infrared types derived here. 
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fast moving solar system objects can appear as detections 
at Y and J while being undetected (at that position) in H 
and K (e.g. LAS YJ observations are sometimes taken on 
different nights to the HK observations of the same region). 
Galaxies can also have unusual colours that could contami- 
nate our selections. 

We used deeper follow-up J- and i?-band photometry to 
identify such contamination, re-measuring the J — H colour 
at a higher signal-to-noise ratio and identifying any extended 
sources. Wo also obtained deeper Y-, K- and z-band pho- 
tometry of a number of our best candidates to provide a 
more complete range of colours with which to assess can- 
didate properties. Wc used a variety of facilities to obtain 
our infrared and optical follow-up photometry, and wc sum- 
marise these measurements for confirmed T dwarfs in Table 
2. In general this table presents photometry for the newly 
discovered T dwarfs, but also includes additional a-band 
measurements for two previously reported LAS T dwarfs 
from Lodieu et al. (2007). 

3.1 Near infrared photometry 

Near infrared follow-up photometry was obtained with three 
facilities. We used WFCAM and the UKIRT Fast-Track Im- 
ager (UFTI; Roche et al. 2003) on UKIRT, both of which 
employ the Mauna Kea Observatories (MKO) J, H and K 
filters (Tokunaga, Simons & Vacca 2002), as well as the 
UKIDSS Y filter (Hewett et al. 2006). We also used the 
Long-slit Infrared Imaging Spectrograph (LIRIS; Manchado 
et al. (1998)) instrument on the William Herschel Telescope 
on La Palma in the Canaries, which uses MKO J and H 
filters, a Ks filter, and a ZLiris filter that is quite similar (al- 
though slightly narrower band; 1.00-1.07/xm) to WFCAM 
Y. It also has a narrow band methane filter {CHaI\ 1.64- 
1.74/L(m) that we used to measure one of the new T dwarfs 
presented here. 

Observations comprised a series of jittered (by a few 
arcsec) image sets, as well as (for WFCAM observations) 
subsets of micro-stepped images that improve the pixel sam- 
pling during processing. Individual exposure times ranged 
from 10s to 216s, and a summary of the different combina- 
tions of micro-step, jitter, repeat, individual exposure time 
and total exposure times, is given in Table 2. All observa^ 
tions were made during photometric conditions with good 
to reasonable seeing (0.6-0.8 arcsec). The data was dark- 
subtracted, fiat fielded and mosaiced using standard ORAC- 
DR routines for the UFTI data, LIRIS-DR routines for the 
LIRIS data, and the Cambridge Astronomical Survey Unit's 
Vista data-flow system (Irwin et al. 2004) for the WFCAM 
observations. Photometry was then performed using typical 
aperture sizes of ~2-arcsec diameter. 

The UFTI observations were calibrated using UKIRT 
Faint Standards, so the photometric system was thus iden- 
tical to the WFCAM system (Lcggctt et al. 2006). The 
WFCAM magnitudes were obtained via flat-file access from 
the WFCAM Science Archive. The LIRIS J-, H- and Ks- 
band observations were calibrated using 2MASS stars as 
secondary standards where J and H magnitudes were con- 
verted onto the MKO system via transforms from Warren 
et al. (2007c). 2MASS uses the same filter as LIRIS, so 
no 2MASS standard K conversion was necessary when de- 
termining LIRIS Ks magnitudes for our candidates. LIRIS 



CH4/ photometry was calibrated by assuming that the aver- 
age H — CH^l colours of 2MASS secondary calibrators was 
zero (see also Kendall et al. 2007). The LIRIS z observations 
were calibrated using observations of AO stars at similar air- 
mass to our target observations, where {z — J)ao=0 was 
assumed. 

In order for the T dwarf photometry to be in the 
same photometric system, we converted our LIRIS z and 
Ks measurements to MKO Y and K. Transforms were de- 
termined via synthetic colours derived from measured T 
dwarf spectra multiplied by the appropriate atmospheric 
and filter transmission profiles, following the methods of 
Hcwctt et al. (2006). Wc derived Y - z^ris and K - 
(LIRIS) colours for the standard T dwarfs from B06 
as well as for the T dwarfs used in Hewett et al. (2006), 
excepting objects that have been shown to be unresolved 
binaries (Liu et al. 2006; Burgasser et al. 2006a). The 
Y — ZLiris colours show no trend with spectral type, and 
a conversion of F=2;Liriji-|-(0.015±0.011) was determined. 
The K — ifs(WHT) colours display a tightly defined se- 
quence monotonically increasing with spectral type such 
that ir=/s's(WHT)-F(0.032xST)-0.01, with a scatter of 0.01 
magnitudes (where ST=1, 2, 3... for Tl, T2, T3...). 

3.2 Optical photometry 

Optical a— band follow-up photometry was obtained with 
the European Southern Observatory (ESO) Multi-Mode In- 
strument (EMMI) on the New Technology Telescope on La 
Silla, Chile. A Bessel 2-band filter was used (EMMI#611; 
zemmi) which has a short wavelength cut-off of 825^im and 
is an open-ended filter, where the CCD sensitivity provides 
the longer wavelength cut-off. The conditions were clear, 
with seeing of '^l arcsec. The EMMI data were reduced 
using standard IRAF routines including a bias and fiat- 
field correction, as well as the removal of fringing effects. 
Multiple images of the same sky areas were then aligned 
and co-added. Aperture photometry was measured using 
^^1.5 arcsec apertures, and wc calibrated this photometry 
using Sloan sources (York et al. 2000) contained in the im- 
ages as secondary standards. This method allowed us to 
derive image zero points with a typical accuracy of ±0.05 
magnitudes. We transformed the Sloan AB photometry of 
our secondary standards into the EMMI system using the 
transformation of Warren et al. (2007a). We then derived 
Zemmi (AB) photometry for the targets, and then trans- 
formed these into the Sloan AB magnitudes assuming that 
a(AB)=ZBMM/(AB)-|-0.2 (Warren et al. 2007a). Our optical 
«— band photometry is thus all in the Sloan AB system. 



3.3 Photometric results 

We have been able to reject 61 of our 103 DRl and DR2 can- 
didates via follow-up photometry, since the J — H colours 
were found to be ^0.4, and the candidates thus likely M or 
L dwarfs. A total of 23 DRl and DR2 T dwarfs have been 
confirmed amongst the remaining candidates, including 11 
of the 15 new T dwarfs reported here (four arc from DR3), 
and 12 previously reported: one by Geballe et al. (2002), one 
by Kendall et al. (2007), one by Warren et al. (2007a), eight 
by Lodieu et al (2007), and one by Chiu et al. (2007). Near 
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Table 2. Follow-up near infrared and optical photometric observations 



Name Instrument Filter UT date ti„t 

& telescope (m=micro-steps, j=jitters, r=repeats)" 



ULAS J0222 


-0024 


UFTI/UKIRT 


J 


2007 Jul 24 


300s 0=5, te^p=60s) 








H 




1800s {j=5, r=6, tei,p=60s) 








K 


2007 Jul 28 


1590s = 5, r = 6, te:rp=53s) 








Y 




600s 0=5, tei:p = 120s) 






EMMI/NTT 


^EMMI 


2007 Nov 17 


2400s (r=4, t<.^p=600s) 


ULAS J0823 


+ 0024 


EMMl/NTT 


zemmi 


2007 Nov 17 


900s (r=2, tea;p='150s) 


ULAS J0837 


-0041 


WFCAM/UKIRT 


J 


2007 May 9 


320s (m=4, j=2, r=4, tea:p=10s) 








H 




800s (m=4, j=5, r=4, texp=10s) 








zemmi 


2008 Jan 30 


2400s (r=4, tea:p=600s) 


ULAS J0859 


+ 1010 


LIRIS/WHT 


H 


2007 Dec 16 


504s 0=9, r=14, te2,p=4s) 






LIRIS/WHT 


CHil 


2007 Dec 16 


450s 0=9, r=5, te:ip=10s) 






EMMl/NTT 


^EMMI 


2007 Nov 17 


600s (r=l, te:ip=600s) 


ULAS J0938 


- 0011 


EMMl/NTT 


^EMMI 


2008 Jan 29 


1800s (r=3, te:cp=600s) 


ULAS J0939 


+ 0016 


EMMI/NTT 


Zemmi 


2008 Jan 31 


900s (r=2, te:ip=450s) 


ULAS J0958 


-0039 


LIRIS/WHT 


J 


2007 Mar 1 


400s (j=5,r=2, t^xp=ms) 








H 




600s 0=5,r=6, tea;p=20s) 








Ks 




800s (j=5,r=8, texp=20s) 








Zemmi 


2008 Jan 30 


2400s (r=4, texp=600s) 


ULAS J1150 


+ 0949 


LIRIS/WHT 


J 


2007 Mar 2 


200s 0=5, te2.p=40s) 








H 




1200s 0=5,r=12, texp='20s) 








Ks 




1800s 0=5,r=19, texv=20s) 








Y 




1200s (j=5,r=6, texp='^Os) 








Zemmi 


2008 Jan 30 


2400s (r=4, tea:p=600s) 


ULAS J1157 


+ 0922'= 










ULAS J1303 


+ 0016 


LIRIS/WHT 


J 


2007 Mar 3 


200s 0=5, texp=40s) 








H 




1200s 0=5,r=12, texp=20s) 








Ks 




1800s (j=5,r=18, texp=20s) 


ULAS J1315 


+ 0826 


LIRIS/WHT 


J 


2007 Mar 2 


200s 0=5, te^p=40s) 








H 




400s (j=5,r=4, te:rp=20s) 








Ks 




1000s 0=5,r=10, tei:p=20s) 








Y 




1000s (j=5,r=5, te2,p=40s) 








Zemmi 


2008 Jan 31 


2400s (r=4, tea:p=600s) 


ULAS J1501 


+ 0822 


LIRIS/WHT 


J 


2007 Mar 2 


200s 0=5, tei:p=40s) 








jj 




400s (i— 5 r— 4 to^r,— 20s') 








Ks 




800s (j=5,r=8, texp=20s) 


ULAS J1505 + 0703 


UFTI/UKIRT 


J 


2007 Jul 23 


540s 0=9, texp=60s) 








H 




3240s 0=9, r=6, texp=60s) 








K 




3240s 0=9, r=6, te2,p=60s) 








Y 




2160s 0=5, r=2, tei:p=216s) 


ULAS J1544 


+ 0819 


UFTI/UKIRT 


J 


2007 Jul 15 


300s 0=5, te^p=60s) 








H 




1800s 0=5, r=6, tei,p=60s) 








K 




1800s 0=5, r=6, te2,p=60s) 








Y 




1080s 0=5, texp=216s) 


ULAS J1547 + 0053 


LIRIS/WHT 


J 


2007 Mar 2 


200s 0=5, teip=40s) 








H 




600s (j=5,r=6, texp=20s) 








Ks 




800s (j=5,r=8, texp=20s) 








Y 




800s (j=5,r=4, texp=40s) 


ULAS J0024 + 0022'' 


EMMI/NTT 


Zemmi 


2007 Nov 16 


2400s (r=4, texp=600s) 


ULAS J0203 


- 0102'' 


EMMI/NTT 


Zemmi 


2007 Nov 17 


1800s (r=3, tea!p=600s) 



"Photometric observations consist of a combination of microsteps, jitters (or dithers), and repeats, of a single 
integration time (texp) that combine to give the total integration time (ti„(). 

A microstep is an offset of a pcirticular number (e.g. N+1/2) of pixels with respect to the current position. 

''T dwarfs from Lodieu et al. (2007). 

°No photometric follow-up was made of this T dwarf. 



infrared and optical photometry for all 15 newly discovered tained to date), and 18 candidates still await follow-up. Fur- 

T dwarfs is shown in Table 3, and their spectroscopic con- ther discussion of the magnitude limited completeness of our 

firmation is discussed in Section 4. One of the remaining follow-up will be presented in Section 6. 
candidates has follow-up photometric colours that suggest 
it is an additional T dwarf (but no spectrum has been ob- 



Table 3. T dwarf photometry and colours. Unless otherwise indicated, near infrared photometry was measured on the MKO system. Uncertainties are indicated in brackets as integer 
multiples of the last decimal place. The last two T dwarfs are from Lodieu et al. (2007) , and here we present additional z-band measurements. 



Name 






Y 


J 


H 


K 


z 




Y -J 


J-H 


H ~K 


z-J 


ULAS 


J0222 


— 0024 


19.87+0.03 


18.71+0.02 


"in no 1 n no 

19.02+0.02 


1 n 1 o 1 n no 

19.18+0.03 


OO cn 1 n r\'7(i 

22. 50+0. U7 


1. 


.16+0.04 


n o 1 1 n no 

-U. 31+0. 03 


n AC 1 n n /I 

-0.16+0.04 


o "70 1 n nT 

3.79+0.07 


ULAS 


J0823 


+ 0024 


19.93+0.15" 


18.57+0.05 


1 o n/^ 1 n 1 oh 

18.95+0.18 


10 CO 1 n oo/) 


OO on 1 n 1 c 

22.80+0.15 


1, 


.od±U.1d 


n on 1 n 1 n 

-0.39+0.19 


n oo 1 n on 


A OO 1 n 1 r* 

4.23+0.16 


TTT AO 


J0837 


— 0041 


1 r» CIA I r\ nob 

iy.o4±(J.Uo 


lo.o2±U.Uy 


1o.dU±U.11 




OO 1 o 1 n r\£i 


1, 


1 o_i_r\ nrt 

.i2±o.oy 


n no i r\ t A 

-0.Uo±U.14 




o c£t 1 n 1 1 
O.DO±U.ll 


TTT AC 


J0859 


+ 1010^ 


iy.uu±u.u7 


i7.oo±U.Uo 


io.Oo±U.Uo 


io.zD±U.iO 




1 


.i2±u.uy 


n ^n 1 n no 

-U.7U±U.Uo 


U.o2±U.lD 




TTT AC 


J0938 


— 0011 


19.83+0.12 


18. 53+0. Oo 


1 r\ nr\ i r» i r\b 

19.0U+0.19 




oo oT 1 n nrr 

22.27+0.05 


1 


on 1 n 1 o 

.30+0.13 


n y( T 1 n on 

-0.47+0.20 




o T yi 1 n no 

3.74+0.08 


ULAS 


.10939 


+ 0016 




i7.9o±().Uo 


TO A -I in"ii6 




01 '7 A 1 n n CT 

21.74+0.05 


1 


.24+0.08 


n A r: i n "i 1 

-U.45±U.ll 




o "vo 1 n nr* 


ULAS 


J0958 


— 0039 




18.95±0.0d 


T n A r\ 1 n in 

iy.4u±u.iu 


"in o 1 n "1 o 1^ 


OO o cr 1 n nn 

22.85+0.09 





.93+0.16 


n ,1 cr 1 n "1 o 

-U.45±U.lz 


n oo \ r\ A T' 

-0.28+0.16 


o nn 1 n "1 n 

3.90+0.10 


ULAS 


J1150 


+ 0949 


19.92±0.08^ 


18.68+0.04 


19.23±0.06 


19.06it0.05^ 


22.44+0.10 


1. 


.24+0.09 


-0.55±0.07 


0.17+0.08 


3.76+0.11 


ULAS 


J1157 


+ 0922 


17.94=b0.03^ 


le-Sl+O-Ol*" 


16.43±0.02^ 


16.24±0.04^ 




1. 


.13+0.03 


0.38±0.02 


0.19+0.04 




ULAS 


J1303 


+ 0016 


20.21±0.17^ 


19.02+0.03 


19.49±0.09 


20.10±0.17'^ 




1 


.19+0.17 


-0.47±0.09 


-0.61+0.19 




ULAS 


J1315 


_1_ (\QOfi 
-\- UoZD 


on nri-i-n ns^ 


lo.oomu.U'^ 


1 Q cn_|_n 1 n 


±y.Duiiiu. ±z 




1 

±, 




-U.D'^IllU. 11 


-u. lumu. ID 


o.yOlllU. 11 


ULAS 


J1501 


+ 0822 


19.70+0.15'' 


18.32+0.02 


18.30+0.06 


18.53+0.10'' 




1 


.38+0.15 


0.02+0.06 


-0.23+0.12 




ULAS 


J1505 


+ 0703 


20.32+0.02 


18.96+0.03 


19.10+0.03 


19.29+0.03 




1 


.36+0.04 


-0.14+0.04 


-0.19+0.04 




ULAS 


J1544 


+ 0819 


19.80+0.05 


18.53+0.03 


18.49+0.03 


18.73+0.03 




1 


.27+0.06 


0.04+0.03 


-0.24+0.04 




ULAS 


J1547 


+ 0053 


19.37+0.06* 


18.32+0.03 


18.45+0.07 


18.21+0.10^* 




1 


.05+0.07 


-0.13+0.08 


0.24+0.12 




ULAS J0024 
ULAS J0203 


+ 0022-'' 
- 0102-'' 










22.77ib0.10" 
22.llib0.06" 










4.61+0.10 
4.07+0.06 



"zemmi converted into zsDSS- 
''Photometry from the UKIDSS LAS. 
°Also has a measured H — CH4I colour of ■ 
'*Ks(Liris) converted into K(MKO). 
"zLiria converted into Y(MKO). 
■''T dwarfs from Lodieu et al. (2007). 



1.05+0.12, typical of mid-late T dwarfs. 
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4 SPECTROSCOPIC OBSERVATIONS 

Spectroscopic follow-up observations were made at five dif- 
ferent facilities, providing a variety of near-infrared spectral 
ranges and sensitivities. Table 4 summarises which facili- 
ties were used to observe the newly confirmed T dwarfs in 
this work (as well as three candidates that were ruled out 
by spectroscopy), and gives details on the wavelength range 
covered, the date of the observations and the total exposure 
times used. A description of these facilities and the reduction 
methods we used is given in the next sections. 



4.1 Gemini/ GNIRS spectroscopy 

The Gemini Near-Infrared Spectrograph (GNIRS; Elias et 
al. 2006) on Gemini South was used to make quick response 
observations, through programme GS-2007A-Q-15. GNIRS 
was used in cross-dispersed mode with the 321 mm~^ grism, 
the 1.0-arcsec slit and the short camera, to obtain 0.9-2.5/im 
R~500 (per resolution element) spectra. The targets were 
nodded three arcsec along the slit in an ABBA pattern 
using individual exposure times of 240s. Calibrations were 
achieved using lamps in the on-telescope calibration unit. AO 
and early F stars were observed as spectroscopic standards, 
either directly before or after the target observations, at an 
airmass that closely matched the mid-point airmass of the 
target, in order to remove the effects of telluric absorption. 
The observing conditions included some patchy cloud, see- 
ing from 0.5-1.0 arcsec, and humidity ranging from 10-50 
per cent. 

Data reduction was initially implemented using tasks 
in the Gemini GNIRS IRAF (Image Reduction and Analy- 
sis Facility) package. Files were prepared and corrected for 
offset bias using NSPREPARE and NVNOISE, and order 
separation achieved with NSCUT. Each order was then me- 
dian stacked at the A and B positions, and a difference im- 
age obtained using GEMARITH. Flat-field correction was 
not necessary since variations across the 6-arcsec slit are 
less than 0.1 per cent, and any flat-field variations in the 
dispersion direction will subsequently be removed when di- 
viding by the spectrum of a standard. S-distortion correction 
and wavelength calibration were performed interactively us- 
ing the telluric star spectra and argon arc lamp spectra, 
with NSAPPWAVE, NSSDIST and NSWAVE. Further re- 
duction was carried out using custom written IDL (Interac- 
tive Data Language) routines. Apertures (1.5-arcsec wide) 
were centred on the spectra at the A and B positions, and the 
sky residuals were fit (and subtracted) using a surface con- 
structed via a series of least-squares linear fits across the slit 
(excluding pixels within the apertures), with one fit for each 
spatial pixel row. Spectra were then extracted by summing 
within the A and B position apertures and combining. The 
target spectra were flux calibrated on a relative scale using 
the telluric standard spectra (after appropriate interpola- 
tion across any hydrogen absorption lines) with an assumed 
black-body function for Teff = 10000 and 7000 K for AO and 
early F tellurics (e.g. Masana et al. 2006), respectively. The 
spectral orders were then trimmed of their noisiest portions, 
and the spectra normalized to unity at 1.27±0.005/im. 



4.2 Gemini/NIRI spectroscopy 

Gemini's Near Infrared Camera and Spectrometer (NIRI; 
Hodapp et al. 2003) was used, on the Gemini North Tele- 
scope on Mauna Kea, Hawaii, through programme GN- 
2007B-Q-26. NIRI was used in the f/6 mode with a 0.75 
arcsec slit, and with the J grism (312.6 1 mm~^) and G0209 
order sorting filter. This produced R~460 spectra over 1.05- 
1.41/im. The targets were nodded 10 arcsec along the slit in 
an ABBA pattern, and calibrations were achieved using ar- 
gon lamp observations (as with GNIRS), and F stars as tel- 
luric standards. The observing conditions were very similar 
to those of the GNIRS observations. 

NIRI data was reduced using tasks in the Gemini NIRI 
IRAF package. NRESID was used to create bad pixel masks 
from the individual flat-field observations. Combined flat- 
fields were created with NSFLAT, and then bad-pixel re- 
moval and flat fielding of the science data was achieved using 
NSREDUCE. Sky removal was then carried out by subtract- 
ing consecutive AB pairs with GEMARITH, and multiple 
spectra combined with NSSTACK. We also minimised any 
pattern noise in our data using a custom written Python 
script, and wavelength calibrated with NSWAVE and cor- 
rected for any S-distortion with NSSDIST and NSTRANS- 
FORM. Spectra were then extracted and calibrated using 
the same methods and software as were used for our GNIRS 
observations. 



4.3 Subaru/IRCS spectroscopy 

We used the Infrared Camera Spectrograph (IRCS; 
KobayEishi et al. 2000) on the Subaru Telescope on Mauna 
Kea, Hawaii. IRCS was used with its camera in the 
52mas/pixel mode, and with the JH grism and a 0. 6-arcsec 
slit. This resulted in R~100 spectra over a wavelength range 
of 1.0 1.6/im. The targets were nodded seven arcsec along 
the slit in an ABBA pattern. Calibrations were achieved us- 
ing an argon lamp, and F2-5 stars were observed as telluric 
standards. The observing conditions were clear with seeing 
of ~0.7 arcsec. 

The IRCS data was reduced following the same reduc- 
tion procedures as were used for our NIRI data, but using 
generic spectral reduction IRAF packages. The only prac- 
tical difference was that S-distortion correction and wave- 
length calibration were done at the same time, using the 
TRANSFORM package. Arc observations were obtained 
several times throughout the night, and we performed wave- 
length calibration with the arc closest in time to each tar- 
get observation. We estimate wavelength residuals (by inter- 
comparing arcs) of ~20A which should not affect signifi- 
cantly our subsequent analysis. 



4.4 UKIRT/UIST spectroscopy 

One of the T dwarfs (ULAS J 1157-^0922) was observed 
at UKIRT using the UKIRT Imager Spectrometer (UIST; 
Ramsay Howat et al. 2004). The HK grism was used with 
the 4-pixel slit, giving a resolution R = 550. Individual expo- 
sure times on target were 240 s, and the target was nodded 
up and down the slit by 12 arcsec. The instrument calibra- 
tion lamps were used to provide accurate flat-fielding and 
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Table 4. Spectroscopic observations. The table contains the new confirmed T dwarfs, as well as (at the bottom) three that we rule out 
with spectroscopy. 



Name 






Instrument 




Wavelength 


UT date 










& telescope 




range 






ULAS 


70999 — 


0024 


NTRT /nomini N 










ULAS 




0024 


NTRT /nprriini N 




1 0^ 1 ZL1 //m 


9 ^on 9007 




ULAS 


J0837 - 


0041 


NlRl/Gcmini-N 




1.05-1.41^m 


5 Oct 2007 


960s (4 X 240s) 


ULAS 


J0859 + 


1010 


NlRl/Gcmini-N 




1.05-1.41^tm 


2 Sep 2007 


960s (4 X 240s) 


ULAS 


J0938 - 


0011 


NIRI/Gemini-N 




1.05-1.41^jm 


2 Sep 2007 


960s (4 X 240s) 


ULAS 


J0939 + 


0016 


NIRI/Gemini-N 




1.05-1.41//m 


2 Sep 2007 


960s (4 X 240s) 


ULAS 


J0958 - 


0039 


GNIRS /Gemini- 


■S 


9-2 hum 


4 Anr 2007 


960s r4x240s^ 


ULAS 


J1150 + 


0949 


GNIRS/Gemini- 


■S 


0.9-2.5//m 


10 Mar 2007 


960s (4 X 240s) 


ULAS 


J1157 + 


0922 


NICS/TNG 




0.9-2.5/im 


9 Apr 2007 


1200s (4 X 300s) 








UIST/UKIRT 




1.4-2. 5Mm 


27 Jun 2007 


5760s (24 X 240s) 


ULAS 


J1303 + 


0016 


GNIRS/Gcmini- 


•S 


0.9-2.5Mm 


10 Mar 2007 


960s (4 X 240s) 


ULAS 


J1315 + 


0826 


GNIRS/Gcmini- 


•S 


0.9-2.5Mm 


4 Apr 2007 


960s (4 X 240s) 


ULAS 


J1501 + 


0822 


IRCS/Subaru 




1-1. 6^m 


1 Jul 2007 


2400s (8 X 300s) 


ULAS 


J1505 + 


0703 


NlRl/Gcmini-N 




1.05-1.41^tm 


22 Aug 2007 


1200s (4 X 300s) 


ULAS 


J1544 + 


0819 


NIRI/Gemini-N 




1.05-1.41/im 


4 Sep 2007 


960s (4 X 240s) 


ULAS 


J1547 + 


0053 


IRCS/Subaru 




1-1.6/im 


1 Jul 2007 


1920s (8 X 240s) 


ULAS J0039 + 


0002 


IRCS/Subaru 




1-1.6/im 


1 Jul 2007 


1920s (8 X 240s) 


ULAS J1301 + 


0023 


IRCS/Subaru 




1-1.6/im 


1 Jul 2007 


1920s (8 X 240s) 


ULAS J2243 - 


0046 


IRCS/Subaru 




1-1.6/um 


1 Jul 2007 


1920s (8 X 240s) 



wavelength calibration. The F5V star HD 96218 was ob- 
served prior to the target to remove the effects of telluric 
absorption, and to provide an approximate flux calibration. 
Relative flux calibration was improved by scaling the spectra 
to the H and K magnitudes measured in the LAS, and fl- 
nally we scaled the spectra at 1.6/ttm to join onto our broader 
wavelength TNG spectra (see next section). 



4.5 TNG/NICS spectroscopy 

We have also obtained a low-resolution (R~50) near-infrared 
(0.9-2.5Mm) spectrum of ULAS J1157-I-0922 with the Near 
Infrared Camera Spectrometer (NICS; Baffa et al. 2001) on 
the Telescope Nazionale Galileo (TNG). The observations 
were made in service mode by the TNG staff, with seeing 
in the 1.2-2.0 arcsec range. NICS is a multi-purpose instru- 
ment equipped with a HgCdTe Hawaii 1024x1024 detector 
with 0.25 arcsec per pixel, yielding a 4.2 by 4.2 arcmin fleld- 
of-view. We employed the Amici mode with a one arcsec slit, 
yielding a wavelength coverage of 0.9-2.5 //m. Four integra- 
tions of 300 sec were taken for the target in an ABBA dither 
pattern. A and B spectra were averaged, differenced, and the 
resulting image flat-fielded. A one-dimensional spectrum was 
then extracted with the IRAF task APSUM. A standard star 
(HIP48971; AO) was observed immediately after the target 
at a similar airmass, and divided through the science target 
to correct for telluric absorption, before flux calibrating with 
a black-body of the same T^s as the AO star. Due to the low 
resolution of the Amici mode, virtually all the Ar/Xe arc 
lines are blended and cannot be easily used for standard re- 
duction procedures. For this reason, wavelength calibration 
was performed in a NICS-speciflc way, using a look-up ta- 
ble (of wavelength against pixel) based on the instrument's 
theoretical dispersion from ray-tracing. This first-pass wave- 
length calibration was then offset to best fit the observed 



spectra of the calibration sources. No attempt was made to 
correct the very modest level of slit curvature. 



4.6 Spectral classification 

Figure 3 shows the spectra of the 15 new T dwarfs. For each 
spectrum we calculated all available near-infrared spectral 
indices from B06; H20-J, CH4-J, H20-H, CH4-H and CH4- 
K. For the GNIRS and NICS spectra we were able to mea- 
sure all five ratios. For the IRCS spectra we were able to 
measure H2O-J, CH4-J and H2O-H. The NIRI spectra al- 
lowed us to measure the H2O-J and CH4-J, and the UIST 
spectrum covered H2O-H, CH4-H and CH4-K. We also com- 
pared the data by eye to template spectra of T dwarfs de- 
fined by B06 as standard T4, T5, T6, T7 and T8 types. The 
results for the new T dwarfs are given in Table 5. Some scat- 
ter is present in the spectral indices given in Table 5 due to 
low signal-to-noise ratio. Consequently, the direct compar- 
ison with templates has been given more weight than the 
spectral indices in the assignment of spectral types, and the 
adopted uncertainty in Table 5 refiects the range in type im- 
plied by this comparison. Note also that in the case of one 
object (ULAS J1150 + 0949), the H2O-H spectral ratio in- 
dicates a significantly earlier spectral type (T3) than all the 
other ratios for this object (T6-7), and we have thus typed it 
as T6.5p. Table 5 also presents distance estimates for the T 
dwarfs based on their spectral type and J-band brightness, 
using the Mj spectral type relation (excluding known and 
possible binaries) of Liu et al. (2006), and allowing for the 
rms scatter in this relation. 

Figure 4 shows the Subaru spectra of the three candi- 
date T dwarfs that we rule out with spectroscopy. Clearly 
none of these objects show the strong H2O and CH4 absorp- 
tion at l.l-1.2yum and 1.3-1.5/Lim that sculpts the J— band 
and ff— band pealcs characteristic of T dwarfs. After allow- 
ing for some residual telluric noise at ~1.4/im, the spectrum 
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"1 1 1 r 

CH4 

I 
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H2O 



H2O 



^ ULAS J1 31 5+0826 (T7.5) 



^ ULAS J1 150+0949 (T6.5) 




ULAS J 1547+0053 (T5.5) 



ULAS J 1303+00 16 (T5.5) 
-mrvr}. .^'^?*>^!^WlM-v 



Ji^f^^ 



ULAS J0939+0016 (T5.5) 




ULAS J0222-0024 (T5.0) 



ULAS J 150 1+0822 (T4.5) 



ULAS J0938-001 1 (T4.5) 



ULAS J 1505+0703 (T4.0) 



ULAS J0823+0024 (T4.0) _ 



ULAS J1 544+081 9 (T3.5) 



ULAS J0837-0041 (T3.0) 




1.5 2.0 
Wovelength (microns) 
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Figure 3. Spectra of the 15 new T dwarfs. The spectra are all normalised to unity at 1.27it0.005/im. Defining features of T dwarf spectra 
are indicated. Also shown are the Y, Yi, J, H and K bands which we use to calculate spectral peaJc ratios (see Section 5). 
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Table 5. Spectral types derived from indices and by comparison with T dwarf templates (following B06). 



Name 


H2O-J 


CH4-J 


H2O-H 


CH4-H 


CH4-K 


Template 


Adopted 


Distance" (pc) 


TTT A Q Tnooo r\r\iA 


U.z (Z \i-0) 


0.368 


(T5.5) 








1 


1 o.UrnU.o 


00- to 


TTT AG TriQOQ _l_ nr\oA 


u.oyu \^o) 


0.696 


(Tl) 












ou- 10 


ULAS J0837 - 0041 


0.419 {T3) 


0.660 


(T2) 








T3 


T3.0±0.5 


54-76 


ULAS J0859 + 1010 


0.058 (T8) 


0.266 


(T7) 








T7 


T7±0.5 


26-36 


ULAS J0938 - 0011 


0.190 (T6) 


0.484 


(T4) 








T4.5 


T4.5±0.5 


53-75 


ULAS J0939 + 0016 


0.288 (T5) 


0.366 


(T5) 








T5-6 


T5.5±0.5 


37-52 


ULAS J0958 - 0039 


0.125 (T6.5) 


0.341 


(T5.5) 


0.391 (T4) 


0.402 (T5) 


0.130 (T6.5) 


T5.5±l 


T5.5±0.5 


58-82 


ULAS J1150 + 0949 


0.087 (T7) 


0.302 


(T6.5) 


0.455 (T3) 


0.230 (T6.5) 


0.032 (>T7) 


T6-7 


T6.5p±0.5 


42-60 


ULAS J1157 + 0922 


0.428 (T3) 


0.544 


(T3) 


0.453 (T3) 


0.772 (T3) 


0.557 (T2) 


T2 


T2.5±0.5 


24-34 


ULAS J1303 + 0016 


0.105 (T7) 


0.324 


(T6) 


0.332 (T5.5) 


0.367 (T5.5) 


0.025 (^T7) 


T5-6 


T5.5±0.5 


60-85 


ULAS J1315 + 0826 


0.034 (T8) 


0.181 


(T8) 


0.227 (T7) 


0.121 (T8) 


0.080 (>T7) 


T8 


T7.5±0.5 


34-48 


ULAS J1501 + 0822 


0.290 (T5) 


0.418 


{T5) 


0.411 (T4) 






T4-5 


T4.5±0.5 


48-68 


ULAS J1505 + 0703 


0.349 (T4) 


0.572 


{T2.5) 








T4±l 


T4.0±0.5 


66-93 


ULAS J1544 + 0819 


0.411 (T3) 


0.520 


(T3.5) 








T3.5 


T3.5±0.5 


55-77 


ULAS .11547 + 0053 


0.217 (T5) 


0.363 


(T5) 


0.331 (T5) 






T5-6 


T5.5±0.5 


44-61 


"^Assumes single objects. 
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Figure 4. Spectra of the three T dwarf candidates ruled out with 
spcctroscoijy. The absorption feature tiiat wc expect to define T 
dwarf spectra are indicated, and are clearly not present in these 
spectra. 

of ULAS J2243-0046 appears to have the form of a blue 
continuum. The other two sources show evidence of H2O ab- 
sorption at 1.3-1.4/Lim, and appear to be M dwarfs. The LAS 
YJ and HK photometry used to select these candidates was 
taken on two separated nights, and variability might explain 
how some M dwarfs can appear to have blue J — H colours. 
ULAS J1301-I-0023 also showed some evidence for possible 
methane absorption (at the ~l-2cr level) via narrow-band 
CH4S- and CH4l-band photometry (Kendall ot al. 2007). Our 
spectrum of this object does not cover the full wavelength 
range of these filters (1.53-1.63/L(m and 1.64-1.74//m respec- 
tively), so we cannot rule out an intrinsically blue methane 



(C HiS-C H4X) colour. However, it may simply be that the 
relatively low signal-to-noise ratio of the narrow-band pho- 
tometry led directly to the measurement of a blue colour. 



5 T DWARF PROPERTIES 

We now consider the overall spectral morphology of the T 
dwarfs using both their spectra and their photometric prop- 
erties. We measured the relative brightness of the Y—, J—, 
H—, and Tf— band spectral peaks by summing flux in several 
bands. In general we used the bands defined by Burgasscr, 
Burrows & Kirkpatrick (2006), to which we also add an ad- 
ditional F— band at a slightly longer wavelength (which we 
refer to as Yi). Our V— , Yi-, J—, H—, and /sT— band spectral 
peaks cover the 1.005-1.045, 1.045 1.085, 1.25 1.29, 1.56- 
1.60, and 2.06-2.10/xm wavelength ranges respectively, and 
are indicated in Figure 5. Table 6 shows the integrated fiux 
peak ratios for the 15 new T dwarfs. We were able to mear 
sure the Yi/J ratio for all our T dwarf spectra, while the 
short wavelength cut-off in the spectral coverage of several 
of the T dwarfs precludes the measurement of the Y/ J ratio. 

To assess how the spectral morphology of the T dwarfs 
may depend on their physical properties, we compare our 
observed properties with the predictions of model atmo- 
spheres. As a primary comparison we use the BT-Settl mod- 
els (which combine with structure models to give the Lyon- 
Settl models shown in Figure 1), generated with version 
15.3 of the general-purpose stellar atmosphere code Phoenix 
(Hauschildt & Baron 1999). These models use a set-up that 
currently gives the best BT-Settl fits to observed spectra 
of M, L and T dwarfs, updating the micro-physics used in 
the GAIA model grid (Kucinskas et al. 2005, 2006). For a 
summary of the important input physics of these models, 
see Warren et al. (2007a). Note that Burgasser, Burrows & 
Kirkpatrick (2006) and Leggett et al. (2007) have also ex- 
amined these effects in the spectral models of Burrows et al. 
(2006; Tucson models) and Marley et al. (in prep.; AMES 
models), and while the different model trends in spectral 
morphology are in broad agreement, we highlight any differ- 
ences between these model sets and the BT-Settl models. 

Figure 5 demonstrates how the BT-Settl model spectra 
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Figure 5. Model BT-Settle T dwarf spectra showing theoretical T^g, logg and [M/H] variations. 
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Table 6. Ratios of spectral peaks. 



Name 




SpT 


Y/J 


Yi/J 


H/J 


K/J 


K/H 


Unusual property? 


ULAS 


70999 — nn94 


T5.0 




0.997 








IN Kjioy JJUool Uiy UliUo \X<XL llifcJ tcLlllCl t_y 


TTT A'^ 


7ns9^ -i- nn9A 


T4 




SQ9 








Noriri3il 


ULAS 


TOK'^i? — nnzLi 

JUOO/ UU'+X 


T3.0 




0.833 








IN \Jl lild/l 


ULAS 


.10859 + 1010 


T7.0 


- 


0.814 


- 


- 


- 


Normal 


ULAS 


J0938 - 0011 


T4.5 




1.195 








Noisy - possibly unusual metallicity 


ULAS 


J0939 + 0016 


T5.5 




0.918 








Normal 


ULAS 


J0958 - 0039 


T5.5 


0.488 


0.708 


0.521 


0.182 


0.349 


Noisy - possibly low Tg£f 
or unusual metallicity 


ULAS 


J1150 + 0949 


T6.5 


0.299 


0.467 


0.428 


0.199 


0.465 


Low-gravity" 


ULAS 


J1157 + 0922 


T2.5 


0.401 


0.467 


0.652 


0.291 


0.446 


Normal 


ULAS 


J1303 + 0016 


T5.5 


0.416 


0.738 


0.449 


0.110 


0.245 


High gravity 


ULAS 


.11315 + 0826 


T7.5 


0.436 


0.723 


0.472 


0.122 


0.258 


Normal 


ULAS 


J1501 + 0822 


T4.5 




0.658 


0.526 






Normal 


ULAS 


.11505 + 0703 


T4.0 




0.710 








Normal 


ULAS 


J1544 + 0819 


T3.5 




0.638 








Normal 


ULAS 


J1547 + 0053 


T5.5 




0.739 


0.437 






High gravity*" 



"Possibly high [M/H] as well. 
'From photometric colour only. 



change with Tctr, gravity, and metallicity. Our baseline prop- 
erties in this figure are Teff=800K, log 5=5.0, [M/H]=0.0, 
and spectral variations for ±100K, Iq'ts ^o's 
in TefT, logg and [M/H] respectively, are shown. The spectra 
have been normalised to unity in the peak of the J-band. 
As has been noted previously for the Tucson and AMES 
models, TefF variations are significant in the wings of the H- 
band peak and in the peak of the K-hwd. The BT-Settl 
models also suggest that the Y-band peak has some TcS 
sensitivity, with increasing and decreasing TefF suppressing 
and enhancing the Y-band flux (with respect to the J— band 
peak) respectively. By comparison the Tucson models also 
exhibit some Y-band variation with T^g, but in the opposite 
sense to that of the BT-Scttl models, and the AMES models 
show little variation with T^s (Leggett et al. 2007). Grav- 
ity has a significant effect on the model iiT-band peak, with 
higher log g causing JT-band suppression and lower log g en- 
hancing the Jf-band flux. Some observational evidence may 
also support this trend, as for example, S Ori 70 (a possi- 
ble T5.5 member of the ~3Myr crOrionis cluster; Zapatero- 
Osorio et al. 2008), and PLZJ93 (a possible T3-5 member of 
the Pleiades open cluster; Casowell et al. 2007) both show a 
strong if— band enhancement, which would presumably re- 
sult from the low surface gravity of members of these solar 
metallicity clusters. 

In addition, the BT-ScttI Y-band flux is suppressed to 
some extent at lower gravity. However, this effect is not pre- 
dicted by the Tucson or AMES models (with both mod- 
els showing an increase in flux for lower gravity), and this 
trend is thus ambiguous when taking all models into con- 
sideration. Metallicity affects the K-haxiA peak in a similar 
way to gravity, where increasing and decreasing metallicity 
causes K-h&nd enhancement and suppression respectively. 
Indeed, K-haxiA fluxes may be quite sensitive to even small 
metallicity differences, as suggested by Liu, Leggett & Chiu 
(2007) when comparing the spectra of HD 3651B (Burgasscr 
2007) and G1570D. Also, metallicity affects the Y-band flux 
peak significantly. This effect is also seen for the other mod- 
els, although it is not clear how the Y-band peak should 



change. The Tucson and AMES models suggest that Y-band 
suppression could occur at higher metallicity, although the 
BT-Settl models suggest that some enhancement may occur 
(at the expense of shorter wavelength flux). 

We do not make direct comparison to the model predic- 
tions, but instead identify T dwarfs whose relative flux peak 
ratios and/or colours are unusual in the context of T dwarf 
properties typical of the solar neighborhood. We compared 
the measurements of the T dwarfs with previous samples 
(from Knapp et al. 2004; Golimowski et al. 2004; Chiu et 
al. 2006), and found that most of the sample (eleven out of 
fifteen) have ratios that appear consistent with the majority 
of solar neighborhood T dwarfs. However, there are some 
notable exceptions. 

• ULAS J0222-0024 and ULAS J0938-0011 appear to 
show Yj enhancement. While our spectral coverage only al- 
lowed us to measure the Y/J flux ratio of ULAS J0222-0024 
and ULAS J0938-0011, the strongly enhanced flux in the Y- 
peak is suggestive of unusual metallicity. However, there is 
no evidence for am enhanced Y flux from the Y — J colours. 
It is possible that the Yi-band (but not the Y— band) could 
be enhanced for these objects. This would be consistent with 
the BT-Scttl model predictions for a high metallicity object 
(see Figure 5), but the Y spectra of these objects are some- 
what noisy and more spectral coverage would be needed to 
properly conflrm this. 

• ULAS J0958-0039 may show some Y-band and H-h&nd 
enhancement, although the somewhat noisy nature of this 
object's spectrum make this result rather tentative. How- 
ever, it is supported by the relatively blue Y — J colour of 
this object. Although note that the Y/J spectral ratio of this 
object is not unusual. Consideration of the models suggests 
that such an effect could arise from lower Tefr or unusual 
[M/H], but it is not currently clear if either or both of these 
effects are responsible. 

• ULAS J1150+0949 shows strong Y and Y suppression 
and also an enhanced if-band flux. The enhanced K/H ratio 
and relatively red H — K colour is suggestive of low gravity. 
The suppressed Y-band flux for this object also points to 
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low gravity when one considers the BT-Settl models. How- 
ever, note that the Tucson and AMES models (which show 
no significant Y-band gravity sensitivity) could invoke an 
increase in metallicity to explain this F-band suppression. 
While it seems simpler to require only one non-typical prop- 
erty to explain the unusual flux peak ratios, it is not possible 
to make firm conclusions at this stage, in light of the model 
ambiguities. ULAS J0859-)-1010 also has suppressed Yj flux 
and unusually red H — K colour, and could have similar 
physical properties to ULAS J1150-I-0949. 

• ULAS J1303-I-0016 shows strong K-hanA suppression 
and very blue H — K and J — K colours, but has a normal Y- 
band flux peak. This is consistent with the predicted model 
trends for a high gravity object. 

• Finally, ULAS J1547-H0053 shows typical Yi /J and H/J 
spectral ratios, but has unusually red J — K and H — K 
colours. Like ULAS J1150-I-0949 this object may have low 
gravity, although it does not show suppressed Yi flux. 

We summarise the above discussion in the last column 
of Table 6, where we indicate possible non-typical properties 
that might explain the observed spectral morphology and 
colours. 



colour of the two T8+ dwarfs is perhaps slightly bluer than 
other mid-late T dwarfs. 

Although only two T8-f dwarfs are currently known, 
their colours suggests that Y — J may be a particularly use- 
ful colour for identification. The J — H colour shows no 
evidence so far for getting bluer beyond T8, however this is 
based on only two objects, and we cannot rule out the possi- 
bility that the small number of T8+ dwarfs discovered so far 
could be unusual. A young (still contracting) brown dwarf 
for instance, would be brighter than an older counterpart of 
the same Tefj , and hence detectable to greater distance. This 
could provide an explanation as to why the first discoveries 
in this new T^s range might have non-typical spectral prop- 
erties. Clearly a larger sample of T8-I- dwarfs is needed to 
properly understand their range of spectral properties. To 
this end, it could prove desirable to extend LAS searches for 
the lowest TcB objects to even bluer Y — J colours, and pos- 
sibly redder J — H. Optical constraints (e.g. non-detections) 
would be particularly important as the near infrared colours 
of candidates become more akin to those of hot stars, but 
the capabilities of UKIDSS and SDSS combine to give great 
potential for such searches. 



5.1 Searching for the lowest Tcff objects 

Figure 6 shows the available colours of our new T dwarfs as 
well as other T dwarfs that have been discovered or recov- 
ered in the LAS (Geballe et al. 2002; Kendall et al. 2007; 
Lodieu et al. 2007; Warren et al. 2007a; Chin et al. 2007), 
including an additional T8-I- object CFBDS0059, recently 
discovered using the Canada-France-Hawaii Telescope (De- 
lorme et al. 2008) and independently identified in the LAS. 
Different T dwarf spectral types are indicated with different 
symbols (see the key in the top right plot). Also shown are 
the sample of DR2 L dwarf candidates (plus signs; see Sec- 
tion 2.1), and the sample of typical brighter sources from 
the LAS (points; see section 2.1) for comparison. Our YJH 
selection box is indicated with a dotted line. 

Examination of the plots shows several interesting fea- 
tures. For the blue T dwarfs (contained in our selection box) 
there is some evidence for a possible trend inY— J colour as 
one moves from the TO-4.5 dwarfs to the T5-8 dwarfs. This 
trend may be continuing into the T8-I- objects, as suggested 
by the rather bluer colour of the T8.5 dwarfs ULAS J0034- 
0052 and CFBDS0059 (filled circles). The J -H colours of 
T dwarfs generally get bluer out to spectral types of T8, 
although the two T8-I- dwarfs do not have very blue J — H. 
The H — K and J — K colours show strong variations irre- 
spective of spectral sub-class, due presumably to variation 
in surface gravity and metallicity, and such variations ap- 
pear to continue into the T8-I- range. As such, these colours 
may not be strong indicators of T8-I- types in themselves, 
although models do indicate that when combined with cer- 
tain spectroscopic information (e.g. Wj index describing the 
width of the J-band peak) such colours (or flux ratios) can 
provide useful T^g constraints (see fig 7 of Warren et al. 2007 
and fig 6 of Delorme et al. 2008). The z — J colour is clearly 
rather scattered for T dwarfs, and there is certainly no evi- 
dence in Figure 6 for T dwarfs becoming redder in z — J at 
lower Tcff . However, it is interesting to note that the z — J 



6 T DWARF NUMBERS 

In this section we discuss the sample of LAS T dwarfs iden- 
tified thus far in DRl and DR2, and place a magnitude lim- 
ited constraint on T dwarf numbers in this 280 sq degs of 
the LAS. We discuss various completeness issues, and make 
approximate corrections to the T dwarf statistics to account 
for these. 

We first consider the completeness with which we have 
followed up candidates from DRl and DR2. Overall we 
have followed up 82 per cent of all our J 5C19.5 candidates 
(see Section 3). Nearly all those that still require follow- 
up are fainter than J=19. Of the 44 initial candidates with 
J ^19.0, 42 have reliable classifications. The follow-up pre- 
sented here combined with the previous observations re- 
ported by KendaU et al. (2007), Lodieu et al. (2007), Warren 
et al (2007a), and Chiu et al. (2007), have resulted in the 
discovery of 23 T dwarfs with J - J/" sJO.l in the LAS DR2 
sky. Of these, 22 T dwarfs have J s;i9.0 (10 with J <18.5 
and 12 with J=18. 5-19.0). Since an extra half a magnitude 
depth should double the surveyed volume, it can be seen 
that our sample shows statistical evidence for a good level 
of completion down to J=19.0. 

The spread in colour of known T dwarfs means that our 
colour selection {J—H ^0.1) will thoroughly probe only cer- 
tain ranges of T spectral class, and for some sub- type ranges 
a fraction of the T dwarfs will be missed. To assess the spec- 
tral type range that our colour selection probes, we used the 
population of optically selected SDSS T dwarf^ (to avoid 
NIR colour selection effects). SDSS T dwarf searches reach 
2=20.4, and should thus be complete across the T range 
{z - J ~3-4) to J ~16.4. There are currently 19 SDSS T 
dwarfs with J <16.4 and J — Hmko >0.1 which have been 
found in 6600 sq degs of sky (see Chiu et al. 2007 and ref- 
erences therein). Amongst these, 14 have spectral type TO- 
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Figure 6. A series of two-colour diagrams sliowing tlie known LAS T dwarfs, including those first presented by Geballe et al. (2002), 
Kendall et al. (2007), Lodieu et al. (2007), Warren et al. (2007a), Chiu et al. (2007), and this work. 



2 and five have spectral typo T2.5-3.5. None have a spec- 
tral type of T4 or greater. Correcting these numbers to a 
LAS J=19 depth and sky coverage reveals that although 
we would not expect to have missed a significant number 
of ^T4 dwarfs, we would expect to have missed ~21 TO-2 
dwarfs and approximately eight T2.5-3.5 dwarfs. Clearly our 
search should be essentially complete for spectral types ^T4, 
but a correction would be necessary if considering spectral 
type ranges that include earlier T dwarfs. Table 7 shows 
how this correction factor changes for different ranges of T 
spectral class, and how the over-all expected numbers of T 
dwarfs would be affected. 

To take advantage of our good level of completeness and 
facilitate simple comparison with the theoretical predictions 
using the methods of Deacon & Hambly (2006) , wc specifi- 
cally consider the spectral type range T4 and later. Allowing 
for the two outstanding candidates we place a constraint of 
19-21 T4-T8.5 dwarfs with measured J-H ^0.1 and J <19 



in the LAS Dr2 sky, and allow for a additional 10% Poisson 
uncertainty on the upper limit to account for the limited 
statistical accuracy of our SDSS sample analysis. 

To account for measurement uncertainties, we have esti- 
mated the number of ^T4 dwarfs with intrinsic J — H ^0.1 
whose colour may have been statistically scattered out of 
the selection region. To do this we estimated (for each ^T4 
dwarf discovered via YJH detection) the probability that 
J~H colours could have been scattered sufficiently rod-ward 
from the LAS measured value, taking into account the 1-a 
uncertainties on their measured colour. Probabilities were 
then summed, and our result suggests that we would only 
expect ~1.5 >T4 dwarfs to be missed in this way. 

In addition, we account for some level of unresolved bi- 
narity amongst the LAS T dwarfs. Unresolved binaries will 
be brighter than a single T dwarf population, and could 
thus be found in a larger volume of sky. A correction is 
necessary if one wishes to compare T dwarf numbers with 



Fifteen new T dwarfs in the LAS 17 



theoretical predictions for the total number of individual T 
dwarfs. Brown dwarf binaries are generally tight systems 
(e.g. <15AU; Reid at al. 2006), of approximately near equal 
meiss components. A volume limited sub-stellar binary frac- 
tion (BF) has been reported with values ranging from 10 
50 per cent via open cluster photometry and high resolu- 
tion imaging (e.g. Burgasser et al. 2003; Pinfield et al. 2003; 
Lodieu et al. 2007), where 



BF ■ 



(1) 



Nb and Ns being the number of binary and single sys- 
tems respectively. If one makes the simplifying assumption 
of binary components of equal brightness (see Burgasser et 
al. 2003 for a more generalised conversion for a range of 
mass-ratio distributions; e.g. Allen 2007), then the relative 
number of binary and single systems in a magnitude limited 
sample will be 



2\f2BF 
1-BF' 



(2) 



and the number of binary systems (compared to the 
total magnitude limited number of systems Nm) will be 



2^2 



2V2 + ^ 



Nrr 



(3) 



For BF= 10-50%, we thus estimate that 24-74 per cent 
of the sources in our magnitude limited sample could be 
unresolved binaries. Of these unresolved binaries systems, 
only ~35 per cent (l/2\/2x 100%) would be included in our 
J ^19 sample if the binary components were resolvable (i.e. 
if the combined binary J-band brightness is >18.25 then the 
individual T dwarf components would have J >19 and would 
be excluded from our sample). However, for the J ^18.25 
unresolved binary systems, we must count both of the T 
dwarf components. Wc thus derive a binarity correction fac- 
tor of 0.76-0.93 (for BF=10-50 per cent respectively) to 
correct our numbers to represent a magnitude limited sam- 
ple of individual T dwarfs. Note that despite a wide range 
in BF, this binary correction factor has a relatively small 
range. 

We also account for spatial incompleteness which can 
result from the way in which optical catalogues are used to 
inform candidate selection. In our initial selection, a nearby 
optical source might result in a candidate being rejected if 
it was close enough to mascarade as an optical counterpart, 
making the candidate appear too bluer in its optical-infrared 
colour. This effectively decreases the search area by some 
fraction that is dependent on the source density and the 
search criteria used. We account for this with a correction 
factor (e) that is the ratio of clear sky (unaffected by opti- 
cal mis-matches) to total sky considered (Atot), and can be 
derived from probability analysis as. 



(4) 



where u is the area excluded by a single optical source 
and N is the number of optical sources in the area consid- 
ered. For 1 dcg^ this reduces to e = (1 — a)" where n is the 
number density of optical sources. Figure 7 shows e plotted 
against source density for our search of the LAS (solid line) 
as well as for previous searches of 2MASS (dashed line, e.g. 
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Figure 7. Predicted area corrections, e = ^'<"^r j^j. excluding 
regions around optical sources as a function oi^ source density. 
The UKIDSS LAS exclusion regions comprise 2 arcsec radii cir- 
cles centred on SDSS sources, and for comparison tiie 2MASS ex- 
clusion regions comprise 10 aresce radii circles centred on USNO 
sources. The source densities inherent in our search and in previ- 
ous 2MASS searches out of the plane are enclosed by open dia- 
monds. 



Burgasser et al. 2004) . The differences result from the differ- 
ent exclusion regions used, which extended out to 10 arcsecs 
from any USNO-A2.0 optical source for 2MASS searches, 
but only extend out to 2 arcsecs from any SDSS source for 
our search of the LAS. The full effect of this incomplete- 
ness depends on the number density of optical sources in 
the regions of sky being considered, and this in turn is de- 
pendent on the galactic latitude. LAS DR2 extends across 
the b =20-80 deg range, although is mostly concentrated in 
b =40-75 deg. The SDSS source number density ranges from 
~20, 000-40, 000 sources per square degree over the LAS DR2 
sky (indicated by open diamonds in Figure 7), and as can 
be seen, this results in a range of e=0.96-0.98. 

For the 2MASS comparison, despite the optical sur- 
veys being shallower than SDSS, the larger exclusion regions 
(with 10 arcsec radii) combined with the greater range of 
galactic latitude searched (reaching jfej — 15) result in the 
e correction having a potentially much greater impact than 
it does for the LAS DR2 T dwarf sample. We chose to fac- 
tor in a 3% increase to our T dwarf numbers to account for 
this spatial incompleteness effect, but note that this small 
correction has no significant effect on our results. 

Finally, wc account for Malmquist bias. We first con- 
sider the effect of an intrinsic scatter in the Mj of T dwarfs 
(classical Malmquist bias). This scatter will cause the J- 
band magnitude limited sample to probe different volumes 
for any fixed spectral type. Those T dwarfs that are intrinsi- 
cally brighter than average (for their type) will be sampled 
from a larger volume than those that are fainter. This will 
bias the mean Mj to brighter values, and effectively result in 
the overall sample coming from a larger volume than if the 
T dwarfs had a non scattered Mj-spectral type relation. As- 
suming a uniform space density and ignoring any effects of 
slope in the substcllar luminosity function (i.e. rM/rfM = 0), 
the fractional increase in the measured immber of T dwarfs 
(0) is given by 



A0 



0.6(7 
loge 



(5) 
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where a is the rms dispersion (assuming a Gaussian 
distribution) of T dwarf absolute magnitudes. We estimated 
this scatter using the sample of T dwarfs with parallax mear- 
surements presented by Liu et al. (2006) combined with their 
polynomial fit (Mj-spcctral type) to the single T dwarfs (ig- 
noring known and possible unresolved binaries). When one 
considers only the subset of T dwarfs with spectral type ^T4 
(i.e. matched to the DR2 sample), the rms scatter from the 
best-fit polynomial is (t=0.29, resulting in an increased vol- 
ume of A(/>/0=16%. We thus decreased our T dwarf numbers 
to account for this bias. 

An additional Malmquist effect results from photomet- 
ric uncertainty near the sample limit, since more T dwarfs 
scatter into the sample from fainter magnitude (i.e. at a 
greater distance, in a larger volume) than scatter out of it. 
However, since we measured accurate follow-up photometry 
for the majority of our initial candidates (down to J=19.5) 
before wo imposed the J=19.0 cut considered hero, the ef- 
fective scatter at the J=19 limit will be small. We conser- 
vatively estimate this scatter to be cr=0.05, and expect a 
resulting Malmquist correction of no more than 0.5% which 
is insignificant. 

Our final magnitude limited T dwarf number con- 
straints take account of the above corrections, and Eissumes 
Poisson uncertainties associated with sample sizes. We thus 
estimate that there are 17±4 single ^T4 dwarfs with J <19 
in the 280 sq degs of Dr2 sky. Table 7 also presents the 
number constraints derived using different lower limits for 
the spectral type range. 

For comparison, we undertook a series of simulations 
based on those presented in Deacon & Hambly (2006). These 
were scaled to take into account the difference in size be- 
tween the full LAS and DR2. The J-band depth limit was set 
to 19, and the F-band depth limit to 20. This F-band depth 
limit is actually 0.2 magnitudes brighter than the 5(t detec- 
tion limit, but this brighter limit effectively takes into ac- 
count the known drop in detection completeness for J=19.7- 
20.2 (Irwin priv. comm.). Unlike the Deacon & Hambly 
(2006) simulations we set no limit on the J!f-band detection, 
since objects not detected in the //"-band will still be identi- 
fied as part of our YJ selection. Also, note that the Deacon 
& Hambly (2006) simulations used a normalisation value 
of 0.0055 objects per cubic parsec in the mass range 0.1- 
O.O9M0 taken from Burgasser (2004). Here we use a normali- 
sation value of 0.0038±0.0013 per cubic parsec calculated by 
Deacon, Nelemans and Hambly (2008). The results of these 
simulations are shown in Table 8, where the simulated late 
T dwarfs were defined as having T^s < 1300if . Mass func- 
tion power-law indices a and x arc also indicated, where the 
simulations assume functions of the form dn/drn ocm^" or 
dn/d(logm) ocm~^ (following Chabrier 2005). 

When comparing the observed number of T4-T8.5 
dwarfs (from Table 7) to the predicted numbers in Tables 
8, and allowing for the associated uncertainties, it is clear 
that it is not possible to place a significant constraint on 
the birthrate. However, this comparison favours a range of 
a between -0.5 and -1.0. Statistically our T dwarf numbers 
are also reasonably consistent (to within 1.5<t) with a = 0.0. 
However, they arc inconsistent with a of 0.5 and 1.0 at the 
2.0 and 2.5a level respectively. This result is reasonably con- 
sistent with the finding of Metchev et al. (2008), who per- 
formed a combined 2MASS and SDSS T dwarf search and 



Table 8. Simulated late T dwarf (defined as those with 
^eff <1300K) numbers for LAS DR2 (see text), for which we as- 
sumed input forms for the initial mass function oi dn/dm oc m~" 
or dn/d{\ogm) oc , and for the birth-rate of h(t) oc e"''*. 



Constant birthrate 


a = -1.0 a; = -2.0 


15±5 


a = —0.5 X = —1.5 


25±9 


a = 0.0 a; =-1.0 


44±15 


Q = 0.5 x= -0.5 


57±19 


= 1.0 X = 0.0 


119±40 


IMF with a = 0.0 or x 


= -1.0 


(3 = -0.1 


54±18 


/3 = 0.0 


44±15 


= 0.1 


33±11 


P = 0.2 


35±12 


P = 1.0 


31±11 



derived a T dwarf space density that was most consistent 
with a = 0.0 (although based on lower number statistics 
than here). However, L dwarf mass function constraints sug- 
gest higher values of a in the 1-1.5 range (Reid et al. 1999; 
Cruz et al. 2007). The sub-stellar mass function of young 
open clusters are generally consistent with a ~ 0.5. How- 
ever, a lognormal function also offers a reasonable fit to 
observation (see the Pleiades mass function in figure 4 of 
Chabrier 2005). The slope of the lognormal mass function 
gets steeper as one decreases brown dwarf mass, and for 
~O.O4M0 the lognormal mass function slope is consistent 
with an a = 0.0 power law. A detailed comparison between 
mass function constraints from L and T dwarfs is beyond 
the scope of this paper, however, it may be that differences 
occur due to the T dwarfs probing a somewhat lower mass 
range than the L dwarfs. This would suggest that a single 
power-law exponent is not optimal when describing the sub- 
stellar mass function in the field. Clearly a greatly improved 
picture of the mass function will emerge from the LAS as 
the survey area grows. 



7 CONCLUSIONS AND FUTURE WORK 

We have discussed our database selection methods and the 
photometric follow-up that we performed to search for late 
T dwarfs and even cooler reff=400-700K objects in the 
UKIDSS LAS. These techniques have allowed us to make 
essentially complete follow-up of our candidates from DRl 
and DR2 down to J=19. We have also followed up some can- 
didates from DR3 and a fraction of our fainter candidates 
from DRl and DR2 (to J=19.5). Using a variety of spectro- 
scopic facilities we have measured the spectra of our best 
candidates, and have spectroscopically confirmed 15 new T 
dwarfs, bringing the total of number of confirmed LAS T 
dwarf discoveries to 28. 

Compared to typical T dwarf properties, one of our new 
T dwarfs may be metal-poor, two may have relatively low 
surface gravity, one may have higher than normal surface 
gravity, and one may be metal rich or have low surface grav- 
ity. These assessments are based on comparisons between 
the spectral morphology and colours of the T dwarfs with 
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Table 7. T dwarf sample statistics. 



Spectral 


No. in our 


Correction to 


Final DR2 


type range 


J < 19 DR2 


include T dwarfs 


corrected 




sky (280 sq 


with J-H >0.1 


number to 




degs) 




J=19 


T0-T8.5 


22-24 


+95±27% 


35±9 


T2.5-T8.5 


22-24 


+40±15% 


24±6 


T4-T8.5 


19-21 


<10% 


17±4 



"Numbers corrected to account for redder ( J — H >0.1) T dwarfs 
outside our selection, LAS photometric uncertainties, unresolved 
binarity for a binary fraction in the 10—50% range, spatial 
incompleteness and Malmquist bias (see text). 



theoretical models trends, and are thus somewhat specula- 
tive. However, all indications are that T dwaxf spectra are 
quite sensitive to their physical properties, and it seems clear 
that a better understanding of these variations would have 
important implications for our ability to model T dwarf at- 
mospheres, and potentially constrain (via spectroscopy) the 
mass, age and composition of T dwarf populations in the 
field (e.g. Pinfield et al. 2006). To this end, it would be very 
beneficial to identify T dwarfs whose properties could bo 
constrained without the need to study their spectra. Both 
T dwarfs and L dwarfs of this type may be found in widely 
separated binary systems (e.g. Gizis et al. 2001; Wilson et al. 
2001; Pinfield et al. 2006; Scholz et al. 2003; Burgasser et al. 
2000; Luhman et al. 2007) for which the primary object can 
be used to constrain system age and composition. G1570D 
and H[D3651B (Burgasser 2000 and Luhman et al. 2007) are 
good examples, being analysed for age and abundance by 
Geballe et al. (2001), Burgasser (2007) and Liu Leggett & 
Chiu (2007). 

The large volume (and large number of T dwarfs) 
probed by the LAS should yield many T dwarfs in such bi- 
nary systems, and in this context we made a basic search for 
possible companions to all confirmed LAS T dwarfs out to a 
separation of ~10,000AU (at the estimated distances in Ta- 
ble 5), by querying the Simbad astronomical database (op- 
erated at CDS). We searched for neighboring objects with 
either a measured spectral type, a parallax, or a high proper 
motion, for which we could obtain distance constraints from 
either parallax or from spectral type and photometry. These 
distance constraints allowed us to rule out (as possible com- 
panions) all the neighboring objects, by comparison with 
the distance constraints of the T dwarfs. However, proper 
motion measurements for all LAS T dwarfs would facilitate 
a more general search for common proper motion compan- 
ions as the number of LAS T dwarfs grows, which could 
include white dwarf companions that can yield useful age 
constraints (e.g. Day- Jones et al. 2008) and are readily iden- 
tified in SDSS (e.g. Kleinman et al. 2004; Gates et al. 2004; 
Eisenstein et al. 2006). 

The number of blue {J-H <0.1) LAS T dwarfs dis- 
covered with our methods shows statistical evidence for a 
good level of completeness down to J=19. This builds on 
our previous work (e.g. Lodieu et al. 2007) for which a good 
level of completion weis only achieved to J ~18.5. The in- 
creased size of the LAS T dwaxf sample allowed us to place 
some statistical constraints on the sub-stellar mass function. 



Indeed, there may be mounting evidence for a steepening 
(decreasing more rapidly with decreasing mass) sub-stellar 
mass function in the field, with best-fit a values of ~ 1 
and between -1.0 and 0.0 for L and T dwarf populations 
respectively. In any event, mass function constraints should 
be greatly improved by larger samples of LAS T dwarfs as 
the survey continues to grow. In addition, constraints on the 
brown dwarf birth-rate should be attainable by building on 
the LAS L and early T dwarf searches mentioned briefly in 
Section 2.1 (e.g. fig 3 of Deacon & Hambly; fig 10 of Bur- 
gasser 2004). For instance, a sample of ~100 L/T dwarfs 
per ATeff=100K range from 1100-1500K (assuming a flat 
birth-rate) would be capable of ruling out an exponential 
(e.g. Tg ~5Gyr) birth-rate at a ~5ct level of significance. An 
appropriate sample of ~400 late L to mid T dwarfs could be 
identified in the full LAS to J=18.5 (for a flat birth-rate). 

The discovery of the T8.5 dwarf ULAS J0034-0052 
(Warren et al. 2007) in the LAS absolutely demonstrates 
that UKIDSS is probing unexplored Tcs regimes beyond 
the previously known late T dwarfs. Having said this, par- 
allax and adaptive optics measurements for this and other 
very late T dwarfs axe important if we are to test for un- 
resolved binarity and better constrain luminosity and Tcff . 
When comparing the colours of ULAS J0034-0052 with all 
available LAS T dwarfs, it is possible that a blue trend in 
Y — J with decreasing Teff may be present. The J — H and 
z—J colours might also change beyond T8, becoming slightly 
redder and bluer respectively, although this needs further 
investigation. These spectral changes would be of great im- 
portance not only to our understanding of cool dwarf atmo- 
spheres, but also to the potential for finding such objects in 
the LAS. If it is the case that roff=400-700K objects become 
bluer in y — J and redder in J — H than T8 dwarfs, it could 
be less problematic to identify such candidates to greater 
J— band depth in multi-band searches of the combined LAS 
and SDSS databases. Indeed, an evolving set of search cri- 
teria takes full advantage of the capabilities offered by large 
general surveys like UKIDSS and SDSS. 
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